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; ABSTRACT 

Context. In the last decade, very extended old stellar clusters with masses in the range from a few 10* to 10* Mq have been found 
' in various types of galaxies in different environments. Objects with masses comparable to normal globular clusters (GCs) are called 

I extended clusters (ECs), while objects with masses in the dwarf galaxy regime are called ultra-compact dwarf galaxies (UCDs). In 

heavily interacting galaxies star clusters tend to form in larger conglomerations called star cluster complexes (CCs). The individual 
star clusters in a CC can merge and form a variety of spheroidal stellar objects. 

Aims. The parametric study aims to analyze how the structural parameters of the final merger objects correlate with the underlying 
CC parameter space. 

Methods. In this work we systematically scan a suitable parameter space for CCs and perform numerical simulations to study their 
(— I . further fate. The varied sizes and masses of the CCs cover a matrix of 5x6 values with CC Plummet radii between 10 - 160 pc and CC 

^ ^ masses between 10^^ - 10** Mq, which are consistent with observed CC parameters. The CCs of the parametric study are on orbits 

I I with galactocentric distances between 20 kpc and 60 kpc. In addition, we studied also the evolution of CCs on a circular orbit at a 

O ■ galactocentric distance of 60 kpc to verify that also extremely extended ECs and UCDs can be explained by our formation scenario. 

Results. All 54 simulations end up with stable merger objects, wherein 26 to 97% of the initial CC mass is bound. The objects show a 
C/3 ' general trend of increasing effective radii with increasing mass. Despite the large range of input Plummet radii of the CCs (10 to 160 

' pc) the effective radii of the merger objects are constrained to values between 10 and 20 pc at the low mass end and to values between 

15 and 55 pc at the high mass end. The structural parameters of the models are comparable to those of the observed ECs and UCDs. 
The results of the circular orbits demonstrate that even very extended objects like the M31 ECs found by Huxor in 2005 and the very 
extended (r^ff > 80 pc), high-mass UCDs can be explained by merged cluster complexes in regions with low gravitational fields at 
' large galactocentric radii. 

, Conclusions. We conclude that the observed ECs and UCDs can be well explained as evolved star cluster complexes. 

OO , Key words. Galaxies: star clusters: general ~ Methods: numerical 
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1. Introduction which have very large radii above 30 pc. These clusters 

^""^ ' were detected by chance as the automatic detection algorithms 

T-H ; Globular clusters (GCs) are very old stellar objects with typical the MegaCam Survey discarded such extended objects as 
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masses between 10^ M© and lO" Mq (corresponding roughly to ]^ background contaminations. Follow-up observations by 
total luminosities between My = -5 to My = -10), having in iMackev et aP dlOOeh . using the ACS camera of the Hubble 



general compact sizes with half-light radii of a few pc. This mor- gp^^^ Telescope (HST), resolved the ECs into stars proving 

H phology makes them easily observable also m external galaxies ^^eir nature as M31 clusters. They also detected a fourth EC 

C?. with modern telescopes (see| Brodie & Stradei1 | 2006i and refer- ^^^^^^ ^31. The M31 ECs have masses of the order of 10^ 
ences therein). Furthe r observations inc reased the number of E Cs 



in 



The Milky Way has a rich GC system containing 150 GCs m31 to 13 (iHuxoret al.ll2()08h . However, iHuxoret all (l20Tl 

([Harris , 1996). Most of them are compact with sizes of a few pc. showed that the previous estimates of the effective radii were 

Only 13 GCs (or 9%) have an effective radius larger than 10 pc. consid erably too large The new size estimate are well below 

Most of these extended clusters (ECs) are fainter than about 30 pc. IChandar et all (12^^41) observed a part of the disks of 

My = -7, only NGC2419, having a half-Ught radius of about the nearby galaxies M81, M83, NGC6946, MlOl, and M51 

20 pc, has a high luminosity of about My = -9.4 mag. Further using HST and found ECs with effective radii larger than 

ECs in the vicinity of the Milky Way have been found in the 10 pc in four of them. M51 showed a veiy high fraction of 

LMC and the Fornax dwarf gal axv (iMackev & Gilmorell2004 eCs in the observed area: 8 of 34 GCs (24%). ECs are now 



LMC and the Fornax dwarf gal axv (Mackev & Gilmordl2004t eCs in the observed area: 8 of 34 GCs (24%). ECs are now 
van den Bergh & MackevI |2004|; [McLaughlin & van der Marej dete cted in all types of galaxies from dwa - fs to ellipticals 
2QQ5I). (e.g. iLarsen & Brodiell200()i iHarris et al.ll200l iLee et al.ll2005l: 



Comp arable objec ts have also been detected around other 
galaxies. iHuxor et aH (l2005h found three ECs around M31, 
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et al.ll2008l: 



Peng et alJl20q6Haues-Santos et alJ[2007t IStonkute 
Georgiev et al."2009^'Da Costa et al."2009'). 

i^Hilker et al. ( 1999) and Drinkwater et (l2000 ll) discovered 
in the Fornax Cluster compact objects with luminosities above 
the brightest known GCs and which were not resolved by 
ground-based observations. These objects have masses between 
a few 10^ Mq and 10^ Mfr, a nd effective radii between reft = 10 
and 100 pc. Drink water et al . (2000) interpreted these objects 
as a new type of galaxy and reflected this i nterpretation in the 
name "ultra-compact dwarf galaxy" (UCD). iBekki et aH (l2001h 
suggested that UCDs are the remnants of dwarf galaxies which 
lost their dark matter halo and all stars except their nucleus. 
Next to the interpretation as a galaxy, UC Ds were also consid - 
ered as high-mass versions of normal GCs (iMieske et al.ll2002h 
or as merged massive complexes of star clusters dKroupall 19981; 
[pellhauer & Kroupa 2002a). 

Many UCDs have been found now. Next to the Fornax 
Cluster, the y have been observed in the galaxy cluster 
Abell 16 89 dMieske et al.l l2004. around M87 in the Virgo 



Cluster ("Hasegan et al.' '2005^ 'Evstigne eva et al.l |2007|) . flie 
Centaurus Cluster (Mieske et al. 2007), the Coma Cluster 
(Madrid et al.l 1201 Ol) . and Abell S0740 (iBlakeslee & Barbed 
I2OO 8). While most known UCDs belong to galaxy clusters, they 
have also been observed in rathe r isolated envir onments, e.g. in 
flie Sombrero galaxy (M104) by Hau et al. (2009D- 

iForbes et'al . (2008) and Mieske et al. (2008) analyzed larger 
samples of UCDs. They find that normal and extended star clus- 
ters and UCDs form a coherent data set where size and mass-to- 
light ratio increase continuously with their total mass and con- 
cluded that UCDs are more likely bright extended clusters than 
naked cores of stripped dwarf galaxies. The marginally enhanced 
mass-to-light ratios of UCDs can b e explained by slightl 
modified initial stellar mass functions dMieske & Kroupall200 
[Dabringhausen et al. 2009). 

High-resolution HST imaging of gas-rich galaxies expe- 
riencing major interactions has resolved very intense star 
formation bursts. The bursting regions are typically located 
within the severely perturbed disks or tidal tails and are con- 
strained to small complexes that contain a few to hundreds of 
young massive star clusters. Exampl es of such systems are the 
knots in the Antennae galaxies ( Whitmore & Schweizerl[T995t 
IWhi tmore et al. 1999), the comp lexes in the NGC7673 star - 
burst (Homeier et al. 2002), M82 (iKonstan topoulos et alj|2009h . 
Arp24 (Cao & Wu 2007), the "b ird's head ga laxy" NGC 6745 
ddeGriis et al. 20 03|), NGC 694 6 dLarsen et al.ll2002h. Steohan's 
Quint et (iGaUagher ;t al.l 1200 ll) . and NGC 922 (iPellerin et al.l 
I2OIOI) . 

The masses of such complexes vary from about 10^ M0 up 
to a few 10** M0. Bastian et al. (2006) observed star cluster com- 
plexes (CCs) in the Antennae with masses of the order !^ 10^ M^ 
and diameters of the order 100 to 200 pc. iPellerin et aP (1201 Ol) 
found young massive CCs with masses between 10^ M0 and 
10^^ Mq and diameters between 600 pc and 1200 pc in the colli- 
sional ring galaxy NGC 922. One of the most extended CCs has 
been observed by iTran et aP (l2003h in the tail of the "Tadpole 
galaxy" UGC 10214. This CC, which has a mass of the order 
10^ Mq, has an effective radius of 160 pc and a diameter of 
ab out 1500 pc. 

iMengel et"an d2008h observed individual young (ailO Myr) 
star clusters associated with CCs in the Antennae and 
NGC 1487. They compared dynamical mass estimates with de- 
rived photometric masses and found them in excellent agree- 
ment, implying that most of them survived the gas removal phase 
and are bound stellar objects. These young clusters are suffi- 
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Fig. 1. Overview of observed efi'ective radii and mass estimates 
of GCs, ECs and UCDs (see Sect.|2]for details). The grey curve 
indicates the general trend of increasing r^ff as parameterized by 
.Dabringhausen et al.. (2008.) . 



ciently s table to be us e d as b uilding blocks for numerical sim- 
ulations. iBastian et aP d2009l) found three 200 to 500 Myr old, 
apparently stable clusters in the Antennae with very high radial 
velocities relative to the galactic disk, indicating that these star 
clusters will most likely become future halo objects. One cluster 
is surrounded by so far unmerged stellar features in its vicinity. 

Since galaxy-galaxy mergers are anticipated to have been 
more common during early cosmological times it is expected 
that star formation in CCs has been a significant star formation 
mode during this epoch. Indeed, the preponderance of clumpy 
galaxies (ElmegreerL2007„ and references therein) indicates that 
early gas-rich galaxies went through an epoch of profuse CC 
formation. We propose a formation scenario for ECs and UCDs 
where massive complexes of star clusters were formed during a 
heavy interaction between gas-rich galaxies, which lead to the 
formation of various kinds of objects via merging of its con- 
stituent star clusters. It has already been shown in previous pa- 
pers that CCs can merge to form a variety of spheroidal stellar 
objects, such as UCDs, ECs, faint fuzzies and possibly dwarf 
spheroidal galaxies (Kroupa 1998; Fellhauer & Kroupa 2002a.^ 
iBekki et aL,2004; .Bilins et al...2009; .Briins & Kroupa 201 1) . hi 
particular, the young UC D W3 is most naturally un derstood to 
be a merged massive CC dFellhauer & Kroupa"2005h. 

In the present paper, we broaden the scope to analyze how 
the structural parameters of the final merger objects correlate 
with the underlying CC parameter space. We systematically scan 
a suitable parameter space for CCs, covering a large range of 
CC masses and sizes. In Section|2]we summarize the status quo 
of available observational data on ECs and UCDs used as input 
for this parametric study. In Section |3] we describe the method 
and the parameters used for the calculations. The results are pre- 
sented in Section|4]and discussed in Section|5] 



2. Observations 

As already indicated in the previous section, GC-like objects 
with eff'ective radii above 10 pc have been found in various envi- 
ronments covering a large mass range. To allow for an overview 
of their parameters, we compiled a list of effective radii and mass 
estimates of GCs, ECs, and UCDs. The respective publications 
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used to compile the data-base are summarized in the following 
subsections. 



2. 1 . GC systems and ECs 

For the Milky Way and its sate llites, we use the catalog of 
resolv ed massive star clusters of McLaughlin & van der Marell 
(l200l . 

The M ilky Way data were complemented by the GC ta- 
bles from IChandar et alj (l2004h for the spiral galaxies M51 
and M81. In addition, t he catalog of GCs i n the lenticu- 
lar galaxy NGC 1380 (Chi es-Santos et alJ l2007h and the faint 
fuzzy star cluster s foun d in the lenticular galaxy NGC 1023 by 
iLarsen & Brodiel (l2000l) were used. The catalog of GCs in 68 
dwarf galaxies dOeorgiev et al. 2009) was added to account for 
GCs and ECs found in dwarf galaxies. 

Next to the aforementione d GC c atalogs miscellaneo us ECs 
fromM31 (Uuxor et al."201 1), M33 (IStohkute e t ar'2008), Scl- 
dEl JPa Costa et al. 2009), and NGC 6822 (Lee et al. 2001) 
we re added to the sample. 

iMcLaughhn & van der Marell (l2005h estimated cluster 
masses for all objects using population-synthesis models to 
define a V-band mass-to-light ratio for every cluster We use a 
typical mass-to-light ratio of M/Ly = 2, which i s consistent 
with the results of iMcLaughlin & van der Marell (l2005l) . to 
estimate masses for the remaining objects of this section. 
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Fig. 2. Three exemplary initial distributions of star clusters 
(small circles with radius R^^) in a CC (surrounding circle with 

radius R'^^^) projected onto the xz-plane. (a) R^^ = 10 pc, (b) 
R^^ - 40 pc, and (c) R^^ - 160 pc. They have a- values of 0.4, 
0.1, and 0.025, respectively. 



2.2. UCDs and massive GCs 



iReikuba et al.l (l2007l) and lTavlor et al.l (l2010l) analyzed a sample 
of massive GCs in the nearby giant elliptical galaxy NGC 5 128 
(Centaurus A). Their sample contains 22 massive GCs with 
masses larger than 10'' M0. 

The parameters of UCD s wer e compiled fo r the 
Fornax Clu ster (Mi eske et all l2008t iRichtler et al. l2005t 
Evstigneeva et al. 2008 ) , the Virgo Cluster (IHa^egan et al.l 
2005:" Evstigneeva et atl l2007i l2008l) . the Centaurus Cluster 



(iMieske et alJ2007l) . and the Coma Cluster (iMadrid et al.ll2010l) . 

In addition, the UCD found in the Sombrero galaxy (Ml 04) by 
[Hau et al. (2009) was considered. 

fBlakeslee & Barbeil (l2008l) identified 15 UCD candidates in 
the ABELL S0740 cluster, which have rather large effective 
radii. We added these UCD candidates to our list, but it should 
be noted that these UCDs are not yet confirmed members of the 
ABELL S0740 cluster and might therefore be background ob- 
jects. 

A large fraction of UCDs has reliable mass estimates. We 
used a typical ma ss-to-light ratio for UCDs of M/Ly = 4 
( IMieske et alj|2008 ) to derive mass estimates for those UCDs, 
where no mass estimate was published so far. 

2.3. Overview ofGC, EC, and UCD parameters 

Figure [1] shows the effective radii of the GCs, ECs, and UCDs 
as a function of their estimated masses. Below masses of lO*" 
M0, the vast majority of clusters have effective radii of a few pc. 
Nevertheless, a few dozens of objects have effective radii larger 
than 10 pc. Most of them have masses of the order 10^ M0. In 
contrast, only very few ECs are found at masses of the order 10^ 

Mq. 

For masses above 10*' M0, there is no clear concentration 
at low effective radii. The radii are more evenly distributed with 
a clear trend of increasing reft with increasing masses. Above 



10^^ Mq all objects have effective radii above 10 pc. The general 
trend of increasing reft, which is added as a grey c urve in Fig. 
[Tj was parameterized bv lDabringhausen et al.l(l2008l) . While this 
line provides a trend for the majority of massive GCs and UCDs, 
the scatter is quite large and a number of objects are located far 
from this line. 

Most UCDs show reft smaller than 35 pc. However, one UCD 
in the Fornax Cluster and one in the Virgo Cluster have effec- 
tive radii larger than 90 pc. These two UCDs show a core-halo 
structure, where the cores have an effective radius of about 10 
pc (Evstigneeva et al. 2007) . In addition, four UCD candidates 
in the Coma cluster (IMadrid et al.l 1201 0|) and the majority of 
the UCD candidates in Abell S0740 ("Bla keslee & Barb'e3l2008l) 
are considerably larger than 35 pc. As most of these objects are 
not yet confirmed as cluster members, they might also be back- 
ground galaxies. 



3. Numerical method and set-up 

3.1. Scenario and model parameters 

The formation scenario described in this paper starts with newly 
born complexes of star clusters covering a large range in masses 
and sizes. We model the dynamical evolution of various CCs 
leading to merger objects. We do not, however, consider the 
galaxy-galaxy interaction, which formed the CCs in the first 
place. 

Individual young massive star clusters were analyzed in de- 
tail by Mengel et al. (2008) and Bastian et al. (2009), leading to 
a combined sample of 25 objects. The median effective radius of 
these 25 young massive star clusters is 4 pc. We use this value 
for our individual star clusters, which are the building blocks of 
the CC models. 

All CC models in this paper consist of A^^ = 32 star clusters. 
The individual s t ar clusters are represented by Plummer spheres 
(|Plummeilll9"TTt iKroupal l2008h with A^,^'= = 100000 particles. 
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Fig. 3. Parameter space for the CC models. The parameters CC 

size, R'"F, and CC mass, M'~''", constitute a matrix of 5x6 values. 
' pi ' ' ' 

All circles mark the 30 simulations on orbit 1 . The open circles 
indicate the 18 additional simulations on orbit 2. The arrows 
indicate the increase of CC-crossing time r^os,, and a- and 
values (see Sect. I3.4l i. The dashed line separates the CC models 
with /? < 1 from the CC models with /3 > 1. CC_34 is an exam- 
ple of the nomenclature of the models. The first index is the row 
number which corresponds to a CC Plummer radius R^^. The 
second index indicates the column synonymous to a CC mass 



The Plummer radius, which corresponds to the effective radius, 
is chosen to be 4 pc for all models. We select a cutoff radius 
of iJf^t = S/Jpf = 20 pc. For each CC, the 32 stai- clusters have 
the same mass, which is 1/32 of the corresponding CC mass. The 
initial velocity distribution of the star clusters is chosen such that 
they are in virial equilibrium. 

The observed CCs sho w a clear concentration of star clus- 
ters towards their centers ( Tran et a"Lll2003t iBastian et al. I l2006t 
[Pellerin et al. 2010). Unfortunately, no detailed observational 
constraints on the distribution of star clusters in the complex and 
their dynamical state are available. In the absence of observed 
density profiles of CCs, we choose a simple model and distribute 
the star clusters in the CC models according to a Plummer dis- 
tribution, which is truncated at the cutoff radius R^^^ - 5R'^I^. 
This cutoff radius is large enough to prevent a clear break or 
edge in the spatial distribution and small enough to avoid single 
star clusters at very large distances that would be stripped away 
immediately. If we would increase our cutoff radius R^^^ from 5 
times the Plummer radius R^^ to infinity, we would have only 
one or two star clusters beyond our actual cutoff radius. Hence, 
the exact value of the cutoff radius will have a negligible impact 
on the results. Figure |2] shows three exemplary initial distribu- 
tions of star clusters in the CC. We used the same seeds for the 
random number generator to generate the same distribution of 
star clusters scaled to the corresponding Plummer radius of the 
CC model. In a previous paper on the e xtended Milky Way clus- 
ter NGC2419 (iBriins & KrouDall20rH) . we have demonstrated 
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Fig. 4. The orbits projected to the xy-, the xz-, and the yz-plane. 
Orbit 1 is a polar orbit and orbit 2 is an inclined orbit between 
galactic radii 20 and 60 kpc. 



that the exact initial distribution of star clusters in an extended 
CC leads to variations in the structural parameters mass and size 
of the order 10 to 20%. The initial velocity distribution of the 
CC models is chosen such that each CC is in virial equilibrium. 
A detailed description of the generation of initial coordinates 
(space an d velocity) for Plum mer models is given in the ap- 
pendix of lAarseth et all (11974 . 

The model parameters constitute a matrix of 5x6 values (Fig. 
with CC Plummer radii of R^^ = 10, 20, 40, 80, 160 pc and 

total CC masses of M^c = 10^ ^ 10\ lO", 10\ 10^■^ lO'* Mq. 
The range of sizes and masses are motivated by the observed 
parameters of ECs and UCDs (see Sect. |2} and they are consis- 
tent with observations of CCs (see Sect. [T]i. For the CC masses 
j^cc _ 2Q6.5 j^jjj iQi.5 additional models with R'^^ - 

240 and 360 pc were considered. 



3.2. Gravitational potential and orbits 

ECs are found near late type disk galaxies, lenticular galaxies, 
elliptical, and dwarf galaxies, while UCDs are predominantly 
found close to giant elliptical galaxies. As the gravitational po- 
tential has a larger impact on the low-mass objects and as we like 
to use the same potential for all computations, we chose an an- 
alytical Milky-Way-like potential consisting of a disk-, a bulge-, 
and a halo component. The coordinate system is chosen such 
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that the disk of the host galaxy lies in the xy-plane. The disk 
is mo deled by a Miyamoto-Nagai potential (iMivamoto & Nagail 
[T97l . 



<l>disk(^,z) = - 



GMd 



(1) 



with Md = 1.0 X 10' ' Mq, fld = 6.5 kpc, and bq = 0.26 kpc. Th e 
bulge is represented by a Hemquist potential (lHemquistlll990l) . 



bulge ~ 



r + a^' 



with Mb = 3.4 X 10'" Mq, and flb 
halo is a logarithmic potential, 



1>halo('-) 



'halo 



(2) 



0.7 kpc. The dark matter 



(3) 



with vo = 186.0 km s"', and rhaio - 12.0 kpc. This set of param- 
eters gives a reasonable Milky-Way-like rotation curve. 

In this paper, we focus on ECs and UCDs located far from 
the galactic disk in the halo of the respective galaxies. As orbital 
parameters for such objects are unknown, we chose a polar orbit 
between galactic radii 20 and 60 kpc for our simulations. These 
values are motivat ed by the projecte d distances to the M31 ECs 
of 13 to 60 kpc (Macke v et al.ll20 06') and the projected distances 
of Fornax UCDs between 8 and 74 kpc (Mieske et al. 2008). 
Figure ID illustrates the chosen orbit, which has an orbital period 
of about 860 Myr 

In our formation scenario, the CCs are most likely formed at 
the peri-galactic passage of the parent galaxy where the impact 
of the interaction is strongest. Therefore, we start our calcula- 
tions at the peri-galactic distance and integrate all models up to 
5 Gyr We stop the integrations at 5 Gyr to save computing time, 
as the structural parameters change only very slightly afterwards. 

To analyze the impact of a polar orbit relative to an inclined 
orbit, we recalculated a subset of our models also on an inclined 
orbit (see Fig.|4] Orbit 2). The orbit is expected to have its largest 
impact on the most extended and lowest mass CC models. The 
additional computations are indicated by open circles in Fig.[3j 
In addition, for the most extended and least massive model 
160 pc and M*^^ = 10^-^ Mq), where the tidal field has 
the largest impact, and for the most extended and most massive 
model (/^p,^ - 160 pc and M'^'^ - IQ^ Mq) complementary cal- 
culations on a circular orbit at a galactocentric distance of 60 kpc 
were performed. 

3.3. Numerical method 

The numerical modeling was performed with the particle-mesh 
code SuPERB0x++ developed by Metz (2008). It is a new C+-i- 
implementation of the FORTRAN particle-mesh code Superbox 
dPellhauer et al.ll2000h using object oriented programming tech- 
niques. SuPERB0x++ makes particular optimal use of modern 
multi-core processor technologies. The code solves the Poisson 
equation on a system of Cartesian grids. 

In order to get good resolution of the star clusters two grids 
with high and medium resolution are focused on each star cluster 
following their trajectories. The individual high resolution grids 
have a size of +80 pc and cover an entire star cluster, whereas the 
medium resolution grid of every star cluster has a size between 
+800 pc and +1200 pc embedding the whole initial CC. The 
local universe is covered by a fixed coarse grid with a size of 



+70 kpc, which contains the orbit of the CC around the center of 
the galaxy. All grids contain 128"' grid cells. 

The galaxy is represented by an analytical potential (see 
Sect. l3.2l i. For each particle in the CC the acceleration from the 
galactic potential is added as an analytical formula to the grid- 
based acceleration computed by solving the Poisson equation. 

3.4. General parameters governing the merging process 

The formation process of the merger object depends on the com- 
pactness of the initial CC. A measure of how densely a CC is 
filled with star clusters for a n equal number Af^*' o f star clusters 
is given by the parameter a dFellhauer et al.ll2002l) . 



■sc 



a - 



(4) 



where R^^ and R^^ are the Plummer radius of a single star clus- 
ter and the Plummer radius of the CC, respectively. In general 
high values of a accelerate the merging process because the star 
clusters already overlap in the center of the CC, whereas low val- 
ues hamper the merging process. Our models cover a-values of 
0.4, 0.2, 0.1, 0.05 and 0.025. High values of a (> 0.1) correspond 
to compact CCs with overlapping star clusters in the center (Fig. 
12^ and b) where the majority of star clusters merge within a few 
Myr (Fig. |5]l. Low values of a (< 0.05) correspond to extended 
CCs (Fig.|2};) where the merging process can take up to several 
hundred Myr The parameter a for the model matrix is shown in 
Fig. |3] increasing from extended to compact CCs. 

Also the tidal field has to be taken into account as it coun- 
teracts the merging process. An estimate of the influence of the 
tidal field on the CC is given by the parameter 
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iCC 



-CC 



(5) 



dFellhauer et al.l20"02h . which is the ratio of the cutoff' radius R^^^ 
of the CC and its tidal radius rp*^. An order of magnitude esti- 
mate of the tidal radius is given bv lKingldl962h . 
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R. 



Mcc 



Mgai,R„(3 + e) 



1/3 



(6) 



where Mqc is the mass of the CC, Mgai.Rp is the galaxy mass 
within Rp, Rp and R^ are the peri- and apo-galactic distances, 
and e - (R^ - /?p)/(/?a + ^p) - 0.5 is the eccentricity of the 
orbit. If the star cluster distribution lies within the tidal radius of 
the CC (J3 < I) the influence of the tidal field on the merging 
process is small. Almost all star clusters will merge and only a 
few are able to escape by chance. However, in case of /? > 1, 
a considerable fraction of star clusters can leave the CC before 
participating in the merging process. The larger the value of /?, 
the larger the impact of the tidal field on the formation process 
of the merger object. The y6-value distribution of the CC models 
is displayed in Fig. [3] The low-mass extended CC models have 
yS- values > 1 . The y6-values increase with CC size and decrease 
with CC mass. 

For a detailed analysi s of merging processes the reader is 
referred to iFellhaueretaD (I2b02i) . 

4. Results 

We carried out 54 different numerical simulations to study the 
influence of varying initial CC parameters. 
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Fig. 5. Time evolution of the merger object of model CC_34 
(M*^*^ - 10^ M0, R^^ - 40 pc). Contour plots on the xz-plane 
displayed at f = 0, 10, 70, 280, 1300 and 5000 Myr The lowest 
contour level corresponds to 5 particles per pixel. The pixel size 
is 5 pc. This yields 0.625 Mq pc"^. The contour levels increase 
further by a factor of 3. 



4.1. Time evolution of the merging process 



To illustrate the merging process, the evolution of model CC_34 
(R'^^ - 40 pc, M*^*^ = 10^ M©) is shown in Fig.|5]as contour 
plots on the xz-plane. The snapshots were taken at t = 0, 10, 70, 
280, 1300 and 5000 Myr The initial CC has an a-value of 0.1. 
At f = 10 Myr (top panel) the merger object is akeady in the 
process of formation. In the course of time more and more star 
clusters are captured by the merger object and the merger ob- 
ject becomes more extended. At f = 70 Myr which is about 10 
CC-crossing times, the vast majority of star clusters has merged 
into the merger object. Another 20 CC-crossing-times later the 
merging process is almost completely terminated and 3 1 out of 
32 star clusters have merged. There is still a close companion 
star cluster as a satellite of the merger object which eventually 
(after f = 1300 Myr) also falls into the merger object. After the 
merging process has been completely terminated the merger ob- 
ject becomes slightly smaller and reaches a stable state within a 
few Gyr. As the structural parameters vary only marginally after 
a couple of Gyr the simulations are terminated at f = 5 Gyr 

The general merging process is very similar for all models, 
but the corresponding time-scale varies considerably. A typical 
timescale for a CC is the crossing time of a star cluster through 
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Fig. 6. (a) Number of merged star clusters versus parameter /3 
for different values of a for orbit 1 after 5 Gyr The horizon- 
tal dashed line marks the number of star clusters A^^ - 32 in 
the initial CC. The vertical dashed line separates the CC mod- 
els with j0 < 1 from the CC models with /? > 1 . (b) Number of 
merged star clusters against CC mass for different CC sizes. The 
horizontal dashed line denotes the number of star clusters in the 
initial CC as in Fig. [6^. 
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The values of the crossing time of the CC cover a range between 
0.3 Myr (CC.16) and 340 Myr (CC.51). The general trend of 
the crossing time is indicated by arrows in Fig. |3] The results of 
our calculations show that on average half of the individual star 
clusters have merged after approximately two-and-a-half 



rCC 
' cross* 



4.2. Number of merged star clusters 

The larger the impact of the tidal field, i.e. the larger the value 
of y6, the smaller is the number of merging events. Fig. shows 
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Fig. 7. Remaining fraction of mass of the initial CC which ended 
up in the merger object vs. the initial CC mass for orbit 1 after 5 
Gyr 
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dius of the initial CC, R'^P, for different CC-masses, M'^^ for 
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the number of merged star clusters as a function of /3. For models 
with P < I, practically all star clusters merge, while for larger 
values of /3 an increasing number of star clusters are able to es- 
cape and align along the orbit. 

Figure|6]3 demonstrates how the number of merged star clus- 
ters depends on the mass and size of the initial CC. The number 
of merged star clusters increases with CC mass and decreases 
with CC size. The number of merged clusters becomes as low as 
13 for the least massive and most extended model (M'^'^ - 10^-^ 
Mq and R^^ - 160 pc). As the unmerged clusters remain com- 
pact GCs, the merging of 13 clusters to one merger objects leads 
to 19 compact GCs originating from the same initial CC. 
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Fig. 9. Effective radii rgff of the merger objects against the 
merger object masses Mend for different CC sizes for orbit 1 af- 
ter 5 Gyr The dashed lines connect models with the same initial 
CC mass. Grey symbols represent additional models. 



4.3. Correlation of structural parameters of merger object 
with CC parameter space 

Naturally the number of merged star clusters has a substantial 
influence on the structural parameters of the merger object. The 
fraction of the merged mass is compared to the initial CC mass 
in Fig.Qfor varying CC sizes. For compact models (a > 0.1), 
where almost every star cluster merges and mass loss is small 
(less than 20%), the final mass of the merger object is compa- 
rable to the initial CC mass. In contrast, the enclosed mass of 
extended CCs strongly depends on the initial CC mass and size. 
For the most extended CC models the merger object masses lie 
between 25 and 80 % of the initial CC mass. The smaller the CC 
mass and the more extended the CC the larger is the influence 
of the tidal field. The CC experiences a larger mass-loss and the 
merging process gets suppressed. 

Figure [8] shows the effective radius r^ff of the merger ob- 
ject vs. the Plummer radius of the CC, R'^^, for different CC- 
masses, M'^'^. The effective radius corresponds to the projected 
half-mass radius, which is the radius within which half of the 
mass is included when projecting the merger object on the sky. 
Compact CCs result in merger objects with effective radii com- 
parable to the Plummer radius of the CC, while extended CCs 
result in merger objects with effective radii that are significantly 
smaller than the corresponding CC Plummer radius. A CC with 
a Plummer radius of R'^^^ = 10 pc leads to merger objects with 
sizes between 10 pc and 15 pc while a CC with a Plummer ra- 
dius of R^^ - 160 pc yields an effective radius range of about 
15 to 55 pc. The more extended the CC becomes the larger is the 
spread in the effective radii of the merger objects. 

For high CC masses of M'^'^ > 10^^ M0 the effective radii 
increase with increasing Plummer radii. For the lower mass CCs, 
the effective radii decrease again for large Plummer radii R'^'^ - 
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Fig. 10. Surface density profiles of the merger objects with R'^^^ Fig. 11. Surface density profiles of the merger objects with 
80 pc and M'^'^ = 10^ \ 10*'■^ and 10^-'* Mq for orbit 1 after 5 = lO^'^ M© and R'^^ = 10, 40, and 160 pc for orbit 1 after 
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4.4. Trends in the r^g vs. Mend space 

The parameter space of the CC models covers the R'^^ vs. M*^*^ 
space uniformly (Fig.O. The corresponding rgff vs. Mend space 
of the merger objects is shown in Fig. |9] For the most compact 
CC models with Plummer radii of R'^^ = 10 pc the effective radii 
and masses of the merger objects are very similar to the Plummer 
radii and masses of the CC, while for the most extended mod- 
els Teff strongly increases with increasing Mend. FigurefTOlshows 
surface density profiles of models with R'^^^ = 80 pc and masses 

of M^'^ = 10^ ^ lO'' ^ and 10^^ Mq. The surface density profiles 
are well represented by King profiles. The structural parameters 
central surface density, core radius, effective radius and tidal ra- 
dius increase significantly with mass. 

For a given CC mass, increasing the size of the CC results 
in a larger mass loss and larger effective radii of the resulting 
merger objects. For the lowest-mass models the effective radii 
decrease again for the largest R'^^^. Figure [TT]illustrates how the 
surface density profiles of the merger objects change with the 
CC Plummer radius (R^^ - 10, 40, and 160 pc) for a CC mass 

of M'"'" = 10^'^ M0. The merger objects show King-like pro- 
files. Increasing R^^ from 10 pc to 40 pc leads to a lower central 
surface density and larger values in the outer parts resulting in 
a larger effective radius. The merger object with R^^ - 160 pc 
suffered a major mass-loss, which leads to considerably lower 
surface densities especially at intermediate radii (5 to 50 pc) re- 
sulting in a lower effective radius. 

For the very extended, low-mass models, the parameter /? is 
much larger than one, i.e. a considerable number of star clus- 
ters of the initial CC is located outside the tidal radius, lead- 
ing to a rapidly decre asing number of merged star clusters. 
iBriins & Kroupal (1201 Ih studied such a turnover in detail for the 
Milky Way EC NGC 2419 and found that the turnover occurs at 
those R^^, where the parameter /3 is sufficiently large to allow 
entire star clusters to escape the merging process. 



As high mass models have larger tidal radii, their size con- 
tinuously increases with CC size up to R^f = 160 pc. However, 
increasing the CC size further will eventually result in decreas- 
ing Teff also for high-mass models. This is demonstrated for 
two additional models for CC masses of M^'^ = 10^-5 M0 and 
j^cc ^ iQi.5 -pjjg QQ si2es were extended to R^^ = 240 

pc and 7?pj^ = 320 pc (grey symbols in Fig. |9]l. For both CC 
masses the results show a clear turnover in the effective radii of 
the merger objects for these large CC sizes. 

The turnover leads to degenerate states in the merger-object 
space, as a relatively compact CC can produce the same merger 
object as a more massive CC having a significantly larger CC 
size. 

The turnover is a general feature of the merging scenario 
which occurs when a significant fraction of star clusters is lo- 
cated beyond the tidal radius of the initial CC. Therefore the 
merging star cluster scenario predicts for each CC mass an up- 
per size limit of the merger objects. The exact CC sizes, where 
the turnover occurs, will also depend considerably on the ini- 
tial configuration, i.e. the exact distribution of star clusters in the 
complex, the number of star clusters constituting the CC, and the 
orbit. 

4.5. Velocity dispersion and dynamical mass 

The global line-of-sight velocity dispersion, cr, is an important 
observable parameter as it can be used in combination with reff to 
es timate the dyn amical mass, Mjyn, of a star cluster. According 
to lSpitzej (Il987h . Mjyn can be estimated by 

Mdyn « 9.75 (8) 

where G is the Gravitational constant. Figure [12] shows the 
global line-of-sight velocity dispersion, cr, as a function of the 
enclosed mass of the merger objects. The turnover, which was 
discussed in the previous section, is clearly seen also in the ve- 
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Fig. 12. Global line-of-sight velocity dispersion cr of the merger 
objects as a function of the enclosed mass Mend for orbit 1 after 
5 Gyr. The solid and the dotted line show the scaling relation 
(Eq. ISJ for objects with an effective radius of 10 pc and 30 pc, 
respectively. 



locity dispersion, as rgfi and cr are not independent. For a given 
mass, an increasing reff results in a decreasing cr. The solid and 
the dotted line in Fig. [12] show the scaling relation of cr versus 
mass for objects with an effective radius of 10 pc and 30 pc, 
respectively. 

For our models, the dynamical masses calculated according 
to Eq. |8] are in very good agreement with the enclosed masses 
with a scatter of about five percent. The small deviations are due 
to slight deviations from virial equilibrium and due to the fact 
that Eq.|8]is only a rough estimate of the dynamical mass. 



4.6. Impact of polar orbit 

We recalculated 18 CC models on an inclined orbit (see Fig. |4] 
Orbit 2) to estimate its impact on the structural parameters of the 
merger objects. In order to save computing time we only recal- 
culated CC models where the inclination of the orbit is expected 
to have a measurable effect. These are the extended CC models 
with large CC-crossing-times, large values of j3, and low values 
of a. These models are indicated by open circles in Fig.|3] 

Figure \T3h compares the effective radii of the inclined orbit 
with those of the polar orbit evaluated after 5 Gyr Both orbits 
produce merger objects with comparable sizes. Only the most 
extended models {R^^ - 160), which are most sensitive to the 
tidal field, show significant deviations between both orbits. 

FigurefTJb compares the enclosed mass of the merger objects 
of the inclined orbit with those of the polar orbit. While there 
are some deviations of Menci between the two orbits for the most 
extended models, the values for both orbits correlate very well. 

The results of the inclined orbit as presented in Fig. [13] 
demonstrate that the inclination of the orbit has no significant in- 
fluence on the overall results and trends of our parametric study. 
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Fig. 13. (a) Effective radius reft,orbit2 of the inclined orbit vs. ef- 
fective radius reff oibiti of the polar orbit for four CC sizes, (b) 
Enclosed mass Menci,orbit2 of the inclined orbit vs. enclosed mass 
Menci.orbiti of the polar orbit for different CC sizes. For the solid 
line the values of orbitl are equal to those of orbit2. 



5. Discussion and conclusions 

We systematically scanned a suitable parameter space for CCs 
and investigated their future evolution. The varied sizes and 
masses of the CCs covered a matrix of 5x6 values with CC 
Plummer radii between 10 - 160 pc and CC masses between 
10^^ - 10^ Mq. The results presented in Sect. |4]demonstrate that 
all simulations end up with stable merger objects, which show a 
general trend of increasing effective radii with increasing mass. 
Despite the large range of input Plummer radii of the CCs the 
effective radii of the merger objects are constrained to values be- 
tween 10 and 20 pc at the low mass end and to values between 
15 and 55 pc at the high mass end. The turnover in the r^s vs. 
Menci space (see Fig. |9]i depends on the mass of the initial CC 
and occurs at larger sizes for higher masses. The turnover leads 
to a degeneracy in the rgfj vs. Mend space of the merger objects. 
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Fig. 14. The diagram shows reij as a function of mass of the ob- 
served ECs and UCDs (circles) and the modeled merger objects 
(stars), combining Figs. \T\ and |9] Two additional models with 
= 160 calculated on a circular orbit at 60 kpc (M'^'^ - 10^^ 

Mq ) are plotted as grey stars. The dotted lines connect 
the additional models on a circular orbit with the corresponding 
models on an eccentric orbit. 



pi 

and 10^° 



i.e. very different CC parameters can result in a comparable final 
merger object. In addition the turnover leads to a higher proba- 
bility for merger objects to have intermediate effective radii. 

Figure [14] shows r^ff as a function of mass of the observed 
ECs and UCDs (circles) and our models (stars), combining 
Figs. [T] and |9] The vast majority of the observed ECs and UCDs 
are located within the parameter space covered by the modeled 
merger objects. Only the very extended objects at Mec ~ 10^ 
M0 and the extremely extended UCDs between Mucd = 10^ 
and 10^ Mq are outside the parameter space covered by this 
study. 

In order to verify that less eccentric orbits would produce 
more extended objects, we calculated two additional models. For 
the most extended models with the lowest and the highest mass 
(/?^c = PC M'^'^ = 10-**-** and 10^ Mq) we calculated the 
evolution on a circular orbit at a galactocentric distance of 60 
kpc. The corresponding merger objects have considerably larger 
effective radii of 29 pc and 82 pc (see Fig. [T4l i than the merger 
objects on the eccentric orbit, which have effective radii of 14 
and 54 pc. Due to the lower gravitational field, the masses of the 
merger objects on the circular orbits are larger than those of the 
eccentric orbits. For the lowest mass model the enclosed mass 
increases from 27% to 50%. These results d emonstrate that very 
extended objects like the M3 1 ECs found bv lHuxor etal] (120051) 
and the very extended, high-mass UCDs can be explained by 
merged cluster complexes in regions with low gravitational fields 
at large galactocentric radii. 

The evolution of ECs in a weak gravitat i onal e nvironment 
has also been studied by iHurlev & MackevI (1201 Cl ). who per- 
formed direct N-body models of extended low-mass star clusters 
incorporating a stellar mass function and stellar evolution. They 
concluded that extended star clusters with an initial mass of 5.8 
10'* M0 are sufficiently stable to survive a Hubble-t ime in a weak 
gravitational field environment. iBekki et al. I (l200l modeled the 
first 70 Myr of the merging of high-mass star clusters without an 



100 
80 

60 
40 



20 

IT 

J 8 
6 



2 
80 

60 
40 



^ 20 



10 

E 



★ 



★ 

9- o o 



★ 
★ 



★ 



-x X 
X , X 



x^^ 



merger objects 
compact UCDs 
extended UCDs 

r.„ = 30 pc 
- scaled r 




10° 



10' 
M [M 



10° 



Fig. 15. (a): Effective radii as a function of mass of the merger 
objects (stars) and of UCDs with observed velocity dispersions 
(crosses for compact UCDs with retj < 10 pc, circles for ex- 
tended UCDs and NGC 2419). (b): Global line-of-sight velocity 
dispersion cr as a function of the mass of the same objects as 
in the above figure. The solid and the dotted line show the scal- 
ing relation for objects with an effective radius of 10 pc and 30 
pc, respectively (see Eq.[8]l. The dashed grey curve indicates the 
general t rend of cr for an increasing r^ff (as parameterized by 
Dabringhausen et al ]l2008i seeFig.D. 



external tidal field. Their finding of triaxial shapes of UCDs is 
most likely related to the relatively short simulation time. Figure 
^illustrates that the merger objects reach a s pherically symmet- 
ric shape only after a few Gyr of evolution. iBekki et al.l (l2004h 
find a general trend of increasing rgfj and velocity dispersions, cr, 
with increasing mass. Their effective radii increase from about 8 
pc at a mass of 4 10^ M0 to values of the order of 20 pc at masses 
of 4 10^ M0. 

Figure [TSl shows only those UCDs and NGC 2419 where 
observations of the effective radius and the global line-of- 
sight velocitj^ dispersion are available (fHasegan etal. 2005t 
Evstigne eva et alJ 2007: Mieske et al. 2008; Bau mgardt et aP 
2009:".Hau et all 12009: Tavlor et al.. .2010) . Figure [T5]i, which 
shows Teff versus mass, demonstrates that most of the very ex- 
tended UCDs with Teff > 20 pc shown in Figure[T4]have so far no 
observed velocity dispersions. Figure [TSb shows observed global 
line-of-sight velocity dispersions of NGC 2419 and UCDs and 
of our merger objects (see Sect. 14. 5t as a function of mass. The 
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models have a steeper relation of cr vs. mass than the observed 
UCDs. This is due to the fact that we modeled solely extended 
objects with > 10 pc, while UCDs show a general trend of 
increasing r^ff with increasing mass (see Fig. [TJ. The mean ef- 
fective radii of the UCDs considered in Figure [T5]increase from 
about 5 pc for masses in the interval lO^" to 10^'^ Mq to about 
16 pc for masses between 10^ ° and 10^^ M0. The grey dashed 
line in Figure [15] which combines Eq.[8]with t he par ameteriza- 
tion of Tetf vs. mass from lDabringhausen et al. is a good 

representatio n of the cr vs. mass relation of the observed UCDs. 
The results o f iBekki et al.l ^004) are closer to the observed val- 
ues as their merger objects have much smaller effective radii than 
our objects. However, it should be kept in mind that a consider- 
able amount of UCDs with large sizes (reft > 20 pc) do not have 
observed velocity dispersions, yet. 

The continuous distribution of CC masses used in our para- 
metric study results in a continuous distribution of masses of 
merger objects. In contrast, the observed masses of ECs and 
UCDs show clear accumulations near masses 10'' Mq and be- 
tween 10^ and 10*^ Mq and a very low number of ECs near 10^ 
M0 (see Fig. fT4l i. A straightforward interpretation on the ba- 
sis of the proposed formation scenario would suggest that the 
mass distribution of ECs and UCDs facilitates direct conclusions 
on the mass spectrum of the CCs, which produced the ECs and 
UCDs. 

An interpretation of the available data on ECs and UCDs 
must be done, however, with great care, as the underlying 
datasets (see Sect.|2]l are highly inhomogeneous and incomplete. 

Due to the limited field of view of the Hubble Space 
Telescope, most extragalactic studies on GCs and ECs cover 
only (a part of) the optical disk of the respective galaxies. The 
ECs discussed in this paper are, however, halo objects located far 
from the optical disk of the galaxies. In the Milky Way, 9 out of 
13 E Cs have galactocentric distances greater than 20 kpc (iHarrisI 
\^6). The only massive EC (Mgc ~ 10'' Mq and reft ~ 20 pc) 
of the Milky Way, NGC2419, is located at a distance of about 
92 kpc. A similar trend has been shown for the other two disk 
galaxies in the Local Group: 12 out of 13 ECs associated with 
M31 and both ECs found in M33 have projected distances well 
outside the optical disks of these gala xies (Huxo r et al. 2008; 
IStonkute et al]l2008t iHuxor et al.ll2009l) . While halo ECs might 
be found by chance in projection to the main body of a galaxy, 
the probability is relatively low: if a survey covers a projected 
area of 20 kpc by 20 kpc and a line-of-sight of +100 kpc is con- 
sidered, the resulting volume, which is covered by the survey, is 
only about two percent of the volume of a sphere of radius 100 
kpc, wherein the ECs would be distributed. 

In addition. iLarsen & Brodi^ (l2000h and lLarsen et al.l(l2002l) 
have discovered a population of ECs co-rotating with the disk 
of the lenticular galaxy NGC1023. These so-called faint fuzzies 
have similar structural parameters as halo ECs and are there- 
fore not easily distinguishable from halo ECs projected onto 
the disk on the basis of imaging data alone. A fair fraction of 
ECs found in extragalactic surveys might therefor e be associated 
with th e disks and not the halos of these galaxies.'Burke rt et al. I 
(l2005h analyzed the kinematics of the faint fuzzies and con- 
cluded that they form a ring-like structure within the galactic 
disk of NGC 1023 and that this ring was probably formed during 
a galaxy-galaxy interaction comparable to the Cartwheel galaxy. 
A detailed discussion of faint fuzzies in the context of merged 
CCs is given in Briins et al. (2009), who demonstrated that the 
observed structural parameters of the faint fuzzies are in excel- 
lent agreement with the merged CC scenario. 



Another reason for incompleteness is the difficulty of dis- 
tinguishing ECs from background galaxies. The GC surveys 
covering 100 galaxies of the Virgo Cluster (Jordan et al. 20()5l) 
and 43 galaxies of the Fornax Cluster dMasters et al.il2010l) ap- 
plied a size limit of retf < 10 pc to reduce the contamination 
of background galaxies . Thereby, t hey e xcluded also all ECs 
from their GC catalogs. Peng et al. (2006) used the same Virgo 
Cluster survey data as Jordan et al. (2005) to search for diffuse 
star cl usters and found e.g . in the galaxy VCC798 ab out 30 ECs, 
where Ijordan et all (l2005b found 21 1 compact GCs. iPeng et aTl 
(l2006h have demonstrated that hundreds of ECs await detection 
in galaxy clusters. Without follow-on spectroscopy to determine 
radial velocities, it cannot be decided whether these ECs are as- 
sociated with the main body or the halo of the galaxies. 

While lower-mass ECs up to Mec ~ lO*" Mq are often 
overlooked in surveys, much effort has been made to detect 
and to analyze stellar objects of M w 10^ Mc?) since the dis- 
cov ery of UCDs in the For nax Cluster by iHilker et all (Il999l) 
and iDrinkwater et al.l (l2000l) . Therefore, the incompleteness for 
UCDs is expected to be considerably lower 

In conclusion, the rapidly increasing number of detected ECs 
and UCDs associated with various types of galaxies in different 
environments offers a new perspective to the process of clus- 
ter formation and galaxy evolution. Since galaxy-galaxy mergers 
are anticipated to have been more common during early cosmo- 
logical times it is expected that star-formation in cluster com- 
plexes has been a significant star-formation mode during this 
epoch. Our work provides a unification of the compact state of 
young star clusters and the extended size of ECs by allowing 
ECs to form from a CC of compact star clusters. If the formation 
scenario suggested in this paper is correct, a statistical analysis 
of the observed ECs und UCDs with respect to the results of the 
models has the potential to shed light on the mass spectrum of 
the initial CCs and thereby on the interaction history of galaxies 
at this cosmologically important epoch. 

Considerably larger, more homogeneous, and more complete 
datasets of ECs and UCDs and more detailed observations of 
CCs are, however, necessary to draw statistically significant con- 
clusions on their origin. 
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